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Cn ■ We discuss the relation between partonic distributions near the phase space boundary 

and Wilson loop expectation values calculated along paths partially lying on the light- 



cone. Due to additional light-cone singularities, multiplicative renormalizability for these 
expectation values is lost. Nevertheless we establish the renormalization group equation for 
the light like Wilson loops and show that it is equivalent to the evolution equation for the 
r~| , physical distributions. By performing a two-loop calculation we verify these properties and 

show that the universal form of the splitting function for large x originates from the cusp 
anomalous dimension of Wilson loops. 
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1. Introduction 

It is well known fT^ that for hard distributions such as deep inelastic structure functions, fragmen- 
tation functions, Drell-Yan pair cross section, jet production etc. near the phase space boundary 
the perturbative expansion involves large corrections which needed to be summed jSISj. This 
region is characterized by the presence of two large scales Q^ and M^ with Q^ ^ M^ S> Aqqj-, 
and at any order in perturbation theory the leading contributions are given by double logarith- 
mic terms (a^ln^^Q^/M^). The leading contributions can be summed and the result is given 
by the exponentiation of the one-loop contribution j^. This exponentiated form holds also after 
the summation of the next-to-leading contributions (a^ln^""^ Q^/M^). The exponent is simply 
modified by two-loop corrections in which the term o?^ In Q"^ /M"^ is absent and the coefficient of 
a^ln^Q^/M^ is proportional to the one-loop beta-function 0. This fact suggests us to apply 
renormalization group approach for summing large perturbative contributions to all orders. 

The universal properties of hard processes near the phase space boundary originate from the 
universality of soft emission, which is responsible of the double logarithmic contributions. Soft 
gluon emission from fast quark line can be treated by eikonal approximation for both incoming 
and outgoing quarks. In the eikonal approximation quarks behave as classical charged particles 
and their interaction with soft gluons can be described by path ordered Wilson lines along their 
classical trajectories [HIIZ]. The corresponding hard distribution is then given by the vacuum 
averaged product of the time-ordered Wilson line, corresponding to the amplitude, and the anti- 
time ordered Wilson line, corresponding to the complex conjugate amplitude. The combination 
of these Wilson lines forms a path ordered exponential W{C) along a closed path C which lies 
partially on the light-cone and depends on the kinematic of the hard process 



W{C) ^ (OlPexp Ugj dz,A^{z) j |0) , A,{z) = A';^{z)X'^ (1.1) 

where A" are the gauge group generator. Here the gluon operators are ordered along the in- 
tegration path C and not according to time. Notice that on different parts of the path C the 
gluon fields are time or anti-time ordered. This expression is different from the usual Wilson 
loop expectation value in which one has ordering along the path for the colour indices of A"^ and 
ordering in time for the gluons fields A'f^lz). In eq. fjl.l|) we have ordering both of colour matrices 
and gauge fields along the path. 

In this paper we study the renormalization properties of W{C) for a path C partially lying 
on the light-cone and show that renormalization group (RG) equation for W{C) corresponds to 
the evolution equation [S] for the parton distribution function or the fragmentation function near 
the phase space boundary. In particular, we discuss the relation between the "cusp anomalous 
dimension" of Wilson loop and the splitting function P{z) for 2; — > 1. 

In sect. 2 we establish the relation between W{C) and the distribution function or fragmen- 
tation function near the phase space boundary. In sect. 3 we perform the two-loop calculation 
in the Feynman gauge and in the MS— regularization scheme of W{C) with the path C fixed 
by the kinematics of hard process. In sect. 4 we deduce the RG equation for W{C) and show 
in sect. 5 that this equation gives rise to the evolution equation for the distribution function 
and fragmentation function near the phase space boundary. Sect. 6 contains some concluding 
remarks. 



2. Wilson loop and structure function in the soft limit 

Consider the deep inelastic process for an incoming hadron of energy E, longitudinal momentum 
Pi and mass M probed by a hard photon with momentum q. In the infinite momentum frame 
we have P+ ^ P_, where the light-cone variables for the hadron momentum P^ = (P+, P^, P_) 
are defined by 



P+ = (P + P,)/V2, P_ = (P - P,)/V2, p^ = (Pi,P2) = 0. 

In this frame xP+ is the "+" component of the momentum of the quark probed by the hard 
photon, where x = — g^/2(P ■ q) is the Bjorken variable. For x ^ 1 the hard photon probes at 
short distances the quark which carries nearly all hadron momentum. In this region the hard 
process is characterized by two scales: the virtuality of the photon Q^ = — g^, which gives the 
short distance scale, and the large distance scale (1 — x)Q'^. 

In this Section we show that, after factorization of leading twist contributions, the partonic 
distributions for x —^ 1 are given in terms of a generalized Wilson loop in p.l|) . To show this we 
first make the leading twist approximation and then we perform the x —>■ 1 limit. 

2.1. Summing collinear gluons 

Due to the factorization theorem jU] the differential cross-section of this process is given in terms 
of the partonic distribution function and the partonic cross section. The partonic distribution 
function is a universal distribution measuring the probability to find in the hadron a parton with 
X fraction of the P+ momentum. To leading twist {Q^ ^ M^) the quark distribution function 
has the following representation ^Oj 

P(a;,/x/M) = J^^ ^^-iP+y-^^p\ ^(^) pexp Ug j' dz^A^{z)\ 7+ ^(0) \P) , (2.1) 

where the vector y^ lies on the light-cone with components 

y^l = {y+,yT,y-) = (o,o,i/_). 

The state |P) describes the hadron with momentum P and mass M. Integration over ?/_ fixes the 
value of the "+" component of the total momentum of the emitted radiation to be (1 — x)P+. 
The corresponding transverse and "— " components are integrated over by fixing y^^ = and 
I/+ = 0. Because of these unrestricted integrations the structure function contains ultraviolet 
(UV) divergences which are subtracted at the reference point fi. The quark field operators ^{y) 
and ^(0) are defined in the Heisenberg representation. The path ordered exponential is evaluated 
along the line between points and y^ on the light-cone and ensures the gauge invariance of the 
nonlocal composite quark-ant iquark operator. 

In what follows we will intensively use properties of Wilson lines which will appear in our 
analysis. We now briefly describe the physical meaning of Wilson line entering into the definition 
(J2.H) . Let us start considering the general form of the structure function of the deep inelastic 
scattering as matrix element of two electromagnetic currents: 

l^,e-^''{P\M^)MO)\P), (2.2) 



Figure 1: One-loop Feyninan diagram contributing to the structure function of deep inelastic 
scattering. The gluon with momentum k is emitted by the quark p and absorbed by the quark 
p' in the final state. We use solid line for quarks, curly lines for gluons, wavy lines for photons 
and dot-dashed line for the unitary cut. 

where 

It is well known that in the leading twist approximation the structure function is defined by 
contributions of Feynman diagrams in which only one of the quarks of the hadron | P) participates 
in the hard scattering ^. One of these diagrams is shown in fig. 1. In the infinite momentum 
frame, in the leading twist limit one can neglect the transverse motion of the quark inside the 
hadron and put the velocity of quark equal to the hadron velocity, 

v^ = PjM. (2.3) 

As we will show in sect. 3, all singularities of the distribution function for x — >■ 1 are absorbed 
by Wilson loop depending on the quark velocity v rather than momentum. 

Consider for example the one-loop contribution of fig. 1 to (j2.2p . In the infinite momentum 
frame the incoming quark has momentum p^ with p_ ^ p+ whereas the recoiling momentum P' 
has component p'_ ^ p'_^_. The vertices for the absorption and emission of the real gluon k are 
given by 

-2pk + iO 2p'k - iO 

respectively. For vanishing quark mass the two vertices contain the following singularities: 

• p— collinear: k+ ~P+, ^- ~ P- 

• p'— collinear: kj^ ~ p'_^, A;_ ~ p'_ 

• soft: k- r^ k+ -^ , 

with krp = 2k^k^. According to the LNK theorem ^T], by summing all diagrams for the structure 
function both the p'— collinear and the soft (or infrared) singularities cancel. We then focus our 
attention only on the p— collinear singularities. There are also diagrams with quark-photon 
vertices "dressed" by hard virtual gluons and quarks having momentum /c+, /;;_, fcx ~ Q- The 
contribution of hard virtual subprocesses can be factorized into the partonic cross section. 

^ There are also contributions initiated by gluons. The corresponding Feyninan diagrams necessary involve real 
emitted quarks whose contribution is suppressed for x ~^ 1. 



Consider the vertex for the absorption of gluon k with gauge potential a^(A;) in the momentum 
representation 

^ 2p'k-iQ ' ^ ' 

where the gluon field carries longitudinal and transverse polarizations. It can be easily shown 
that when the momentum k is collinear to p, the components of afj_{k) transverse to k contribute 
to higher twist and can be neglected. Thus, only longitudinal polarization survives and the gauge 
potential becomes pure gauge field 



yk — iO 



This property allows us to sum to all orders the contribution of collinear gluons (leading twist 
approximation) using Ward identities. Inserting this gauge field in (12. 4j] we obtain that the gluon 
emission is factorized into 

^'■^/(0^(^^)=^'''^/'^'''^^^''^'^'^''^=^'''^/°°'^^^^'^'^^^^' ^^-^^ 

where 

poo 

Jf,{x) = y^, / dr 6^{x - yr) 
Jo 

is the classical non-abelian eikonal current of the scattered quark p' in which we neglect higher 

twists contributions of order p'_^_/p'_. This factorization property implies that the quark field 

operator ^(0) describing the scattered quark p' in ()2.2p can be written in one-loop as 

fr(O) = fro(o) fl + tg j^ dz^A^{z) + 0{g'')\ (2.6) 

with integration path along 7/_ axis. Here, 'I'o(O) is a free quark operator and A^{x) describes 
p— collinear gluon. Note that this relation is valid only within the matrix element of currents in 
fl2.2|) and not as an operator identity. 

The generalization of the one loop result ()2.6|) is based on the following physical arguments. 
After the hard scattering the quark with momentum p' starts to emit gluons and quarks moving 
in the same direction. All these particles form p'— collinear jet propagating through the cloud 
of p— collinear gluons created by the incoming quark. Since p— collinear gluons are pure gauge 
fields they cannot change the state of the scattering quark but only modify its phase. Such a 
pure phase factor is given by Wilson line along the light-cone direction y defined in fj2.1|) . The 
generalization of ()2.6|) is then 

^(0) = fp._coi(0)<l'j,[0,oo;A], (2.7) 

where ^y{zi, Z2) is the path ordered exponential calculated from point zi to Z2 along the direction 



y 



$y[zi, Z2] A] = Pexp i^ / dz^Af'i^z 



22 



In (12. 7j) the phase $y[0, 00; A\ describes the interaction of p— collinear gluons with the scattered 
quark while operator f'p/„coi(0) describes the jet of p'— collinear particles produced by the scat- 
tered quark. 



For the quark field operator "^{z), which describes absorption of p— collinear gluons by quark 
with momentum p' in the final state, one obtains 

^{z) = ^-y[oo, y; A] ^p>-coiiz) , (2.8) 

with y = (0,0, ?/_) and ?/_ = z_. After substitution of (|2.7p and (|2.8|) into (j2.2|) we get that the 
full phase of the scattered quark is given by 



<t>^y[oo,y;A]%[0,oo]A] = V exp lig / dz^Af^iz] 



(2.9) 



Moreover, the contributions of p— collinear and soft particles to the structure function can be 
factorized into the distribution function 1)2.111 . The contribution of p'— coUinear jet is factorized as 
usual into the partonic cross section ^. Thus, the Wilson line in the definition of the distribution 
function ()2.H) is the result of the interaction of scattered quark p' with the gluons collinear to 
incoming quark p. Here, we have used the following properties of the Wilson lines: 

• hermiticity: $^[0, oo; A] = $_,y[oo, 0; A] 

• causality: ^y[b,c; A]^y[a,b; A] = ^y[a,c; A] (2.10) 

• unitarity: ^l[a,b; A]^y[a,b; A] = 1 

Notice, that the approximation (|2.9|) is valid in the leading twist limit for an arbitrary values of 
the scaling variable x. 

2.2. Soft approximation for x ^ 1 

For X -^ 1 the total "+" momentum of the emitted radiation (1 — x)P+ vanishes and the 
expression ()2.H) can be further simplified. In the leading twist order only p— collinear and soft 
gluons contribute to ()2.1|) . As already recalled, due to the LNK theorem p— collinear singularities 
survive in parton distribution ()2.1|) . Although soft divergences cancel, their finite contribution 
completely determines the asymptotic parton distribution function for x — *> 1. To show this we 
observe that for a; — > 1, the "+" component of the momentum of any emitted real gluon q is 
positive and vanishing as g-|_ ~ (1 — x)P+. Therefore, p— collinear gluons are only virtual since 
their momentum q has component g+ ~ P+. It means that real emitted gluons are only soft 
(with momenta g+, g_, gx ~ (1 — 3;)P+) while virtual gluons can be either soft or p— collinear. 

These conditions imply that before hard scattering by photon probe, the incoming quark with 
momentum p decays into a jet of p— collinear virtual gluons and quarks emitting and absorbing 
soft gluons. The contribution of real and virtual soft gluons can be factorized from the distribution 
function as follows. We note that soft gluons interact with the p— collinear quarks and gluons, 
forming the p— collinear jet, via eikonal vertices. In this approximation the incoming quark 
with momentum p behaves as a classical charged particle moving with velocity f^ defined in 
()2.3|) and all effects of interaction of the p— collinear jet produced by incoming quark with soft 
gluons are absorbed into phase factors similar to path-ordered exponentials found in the collinear 

^As notices in refs. E!i the latter contribution become large for x — > 1. That is why performing the 
factorization of the structure function for a; ^- 1 one should treat the p'— collinear jet separately from the partonic 
cross section. 



approximation. Therefore, in the x —>■ 1 hmit one can replace quark field operators in the 
definition p.lj) as 

^{y) = ^p.,oi{y)<l>~v[y,oo-A], $(0) = <l>4oo,0;A]^p„coi(0), (2.11) 

where the phase factors are evaluated along a trajectory of a massive classical particle with 
velocity v. Here, the gauge field operators Af^{x) describe soft gluons and the quark fields 
^p-coi{y) and \E'p-coi(0) describe the incoming quark in the initial and final states dressed by 
p— coUinear virtual corrections. 

We did not consider till now the contribution of quark emissions to the distribution function. 
Repeating the previous arguments one finds that for x — *> 1 real quarks are soft. But it is well 
known from power counting |E] that soft quarks give vanishing contribution to the leading twist 
order. Therefore all emitted quarks must be virtual and they contribute either to p— collinear 
jet or renormalize soft gluon self-interaction vertices. These last contributions will be explicitly 
evaluated in subsect. 3.2. 

It should be noticed that similar phase factors ()2.7|1 and ()2.8|1 describe interaction of scattered 
quark with p— collinear gluons and incoming quark with soft gluons, respectively, although the 
two subprocesses are different. This is due to the fact that in both cases gluons do not change the 
quark state: the coUinear gluons are pure gauge fields with vanishing strength; the soft gluons 
interact with quarks via eikonal vertices which preserve the momentum and polarization of the 
quarks. 

Combining the phase factors ()2.9|1 and ()2.1H) from collinear and soft approximations we obtain 
the following representation for the distribution function in the x — > 1 limit: 

F{x,ij/M) = H{ij/M) [ " P^'^e'P+y-i^-)w{Cs) (2.12) 

with the function H{fj,/M) taking into account contributions of p— collinear quarks and gluons 
and 

W{Cs) = {0\<^^,[y,oc;A]<l>y[0,y;A]<^4oo,0;A]\0) = (0|Pexp Ug I dz^A^^iz)) |0) , (2.13) 

where the index S indicates that the photon probe has space-like momentum. The integration 
path Cs = £iU£2U£3 is shown in fig. 2. The rays ii and £3 correspond to two classical trajectories 
of a massive particle going from infinity to with velocity v and from y to infinity with velocity 
—V, respectively. The segment £2 = [0,y] lies on the light-cone. As explained before, since we are 
dealing with a cross-section rather than with an amplitude, the gauge fields are ordered along 
the path rather than the time. Path- and time-orderings are related as follows 

P = T for £1 , P = r for £3 , P = T for £2 and y_ > , V = T for £2 and V- < , 

where T and T denote time- and anti-time ordering, respectively. This path has two cusps at 
the points and y where the direction of the particle is changed by the hard photon probe. 

Thus the quark distribution function for x ^ 1 is given by the product of two functions: 
the Fourier transform of the Wilson loop expectation value, which takes into account all soft 
emissions and gives the asymptotic distribution function for x — >■ 1, and times the function 
H[ix/M) getting contribution from virtual p— collinear quarks and gluons. This is the reason 
why H{fi/M) does not depend on P+(l — x) or, equivalently, on y^. Although F{x,fi/M) cannot 
be calculated in perturbation theory we will find its evolution with the renormalization point fi 
using renormalization properties of W{Cs)- 




Figure 2: Integration path C5 = £i U £2 U £3 for the Wilson loop W{Cs) corresponding to the 
structure function for large x. The ray ii is along the time-like vector f^ from — cxo to 0; the 
segment £2 is from point to y along the light-cone; the ray £3 is from the point y to — 00 along 
the vector ~v^. This path has two cusps at points and y where the quark undergoes hard 
scattering. 

2.3. Fragmentation function in the x ^ 1 limit 

The Wilson loop representation for the structure function in ()2.12|) can be easily generalized for 
the fragmentation function D{x, fi/M) in the limit x ^ 1. This function gives the probability of 
finding a hadron in a quark jet. To the leading twist approximation the fragmentation function 
is given by [TU] 

^g-.p+._/. J2^p^ /|^(y)$,[oo, y; A] |0)7+(0|<l>_,[0, 00; A]^{0)\P, f) 

-00 r 

(2.14) 
where P^ = (P+,0,P_) is momentum of the observed hadron, y = (0,0, y_) is a vector in 
coordinate space which lies on the light-cone, and / denotes the associated final state radiation. 
In the infinite momentum frame with P+ ^ P_ integration over y_ fixes the "-I-" component of 
the total associated radiation momentum to be [l/x — 1)P+. As in Subsection 2.1, the phase 
factors in eq. ()2.14j) are obtained from the resummation of gluons interacting with the recoiling 
quark P' and collinear to the observed hadron P. 

In the limit x — *> 1 the associated radiation becomes soft and we can perform the eikonal 
approximation as in the previous section. Namely, the dependence of gauge fields in the quark 
and anti-quark field operators in ()2.14|) factorizes into two phase factors as in eq. ()2.11|) . In the 
soft {x —>■ 1) limit we obtain the following representation for the fragmentation function 

D{x,ij/M)=H{i2/M) [ °^ P^'^e^P+y-^^-y-)w{CT). 

J-00 27r 

The generalized vacuum averaged Wilson loop operator W{Ct) is given by 

W{Ct) = J2{m.[y,oc;A]%[oc,y;A]\0){0\<^_y[0,oo;A]<l>_,[cx,,0;A]\f) 
f 
= {0\^4y, 00; A]%[0, y; A]$_„[oo, 0; A]\0) 



{0\Vexp[tg(b dz^A^iz)] |0) 

'Ct 
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y 

Figure 3: Integration path Ct corresponding to the fragmentation function for large x. 

where we have used completeness condition for the associated radiation 1 = Xlf l/)(/l ^^"^ 
causality (j2.1Up for the phase factors. The index T indicates that the photon probe has time-like 
momentum. 

The integration path Ct = £i U £2 U £3 is shown in fig. 3 and three parts £j have meanings 
similar to that described before for the space-like process. In this case path- and time-orderings 
are related as follows 

P = T for £1 , P = T for £3 , P = T for £2 and ?/_ > , V = f for £2 and i/_ < 

This path has two cusps at the points and y where the direction of the particle is changed by 
the hard photon probe. 

2.4. One-loop calculation 

We perform the one- loop calculation of W{Cs) in order to explain in this formulation the nature 
and origin of double logarithms and to discuss the analytical properties of W{Cs) in the ?/_ 
variable. 

We parameterize the integration path C^ = ^1 U £2 U £3 = {z^{t)] t G (— C)0, +00)} as follows 



z,{t) 



The relevant diagrams are given in fig. 4 plus the symmetric ones. We perform the calculation 
in the Feynman gauge and in coordinate representation by using the dimensional regularization. 
As we shall see in the sect. 3 this representation is very convenient for the calculation of W{Cs) 
to two loops. Feynman rules in coordinate representation are summarized in the Appendix A. 

The path from point to y consists actually of the product of paths from to 00 and from 
00 to y as shown in fig. 5 for the diagram of fig. 4a. The gluon is either real (fig. 5a) or virtual 
(fig. 5b). In both cases the gluon is emitted at point zi = vti with ti < 0. For real gluon we have 
Z2 = yt2 with 1 < ^2 < 00, while for virtual gluon we have Z2 = yt2 with < 1(2 < 00. The two 
contributions are given by 

dt, / dh D'^''{Z2 -z,)+ dt2 Dl''{z2 - z,) \ (2.15) 

00 Uo Joo ) 



y^^t, 


-00 < i < 
< t < 1 




y^,-Vf,{t-l), 


1 < t < 00 






(a) 



(b) 



(c) 



Figure 4: One- loop diagrams contributing to W{Cs)- Here, the double line represents the in- 
tegration path in the Minkowski space as in fig. 2. The Feynman rules for these diagrams are 
given in the Appendix. 





Figure 5: Diagrams corresponding to fig. 4a for virtual and real gluon (see ()2.15|) ). Due to a 
partial cancellation of Wilson lines, the sum of these two diagram gives the single contribution 
of fig. 4a. Similar cancellations hold for all diagrams. 



where Cp is the quadratic Casimir operator in the quark representation. For y_ positive the 
is time-like {z'^ > and zq > 0) so that cutted and full propagators coincide: 



vector z = {z2 — zi 
Dl^'iz) = D^\z) = 



—g^^D[z). Then, the sum of the two diagrams of fig. 5 is 
W^i^) = g'^Cpivy) f dti [ dh D{yt2 - nh) 

J-oo Jo 



(2.16) 



The same result is obtained for ?/_ negative by deforming the integration path ii. To show this 
observe first that the singularities of the integrand in the ti— complex plane lie in the lower half 
plane. Therefore we can deform the integration path from — oo<ti <OtoO<ti <cxd. After 
this deformation the vector yt2 — vti becomes time-like and we have the same situation as for y_ 
positive. 

Both the real and virtual contributions in ()2.15|) have infrared (IR) and UV singularities. 
The real and virtual IR singularities, coming from ^2 -^ oo, cancel in ()2.16|) . giving a bonded 
range for t2- This cancellation between the two paths from point y to oo is general since it is the 
result of causality of Wilson lines in (j2.1Up . For this reason in fig. 4 we neglect these parts of the 
integration paths. 

Integral in eq. ()2.16|) is ultraviolet divergent and we use dimensional regularization. Using the 
D— dimensional propagator given in the Appendix we obtain 



wP 



9^ 

47r«/2 



dti / dt2{-tl + 2{vy)tit2 + t0y' 

oo Jo 



D/2 
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which is divergent for D = A. To see the origin of these divergences we change the integration 
variable to ti = —2{vy)t2t[ and obtain 

O'' ' Jo Jo 

For D ^ 4 we find two singularities: for ^2 ^ and t'^ ^ 0. In the first case we have Z2 — zi — > 0, 
which corresponds to the integration near the cusp at point 0. This singularity is called "cusp 
singularity" [12]. For t[ ^ we have ^i — >• and {z2 — Zi)^ — ^ 2;| = which corresponds to the 
light-cone (coUinear) singularity. The final unrenormalized result for this diagram is 

^a. . -^CH«.(.. . . - .or^ n3-gmz.-3) ^^^^, 

For the diagram of fig. 4b we obtain 

/O poo 

dh / dt2D4y-v{t2 + h)) 
-oo Jo 

= -i^/^'"''rp/2-l)C7^ / dt^^dt2 
X [2{y_ - v.{t2 + ti) - iQ){v+{t2 + ti) + i^f-"'^ 
This integral for D ^ 4 has an infrared singularity for Z2 — Zi -^ oo. Taking D > A one gets 

,.« . ^CH2>.,« . . - .or. r(3-W(z>-3) ,^^^, 

For the diagram of fig. 4c we have 

^i'^ = -^li'-^T{D/2-l)CF j^ dhj dt2{-{h-t2)' + t0y-''/' 

= -i|^^H/^^^^^)^-^/^^^-^^^^^{^/^rft,t?-^. (2.19) 

In this case we have an IR singularity for D < 4 as ti — ^2 — >^ oo and an UV divergence for D > 4 
as ti — *> and ti — ^2 — ^ 0. The IR divergence of Wc is cancelled by the IR divergence of W^ . 
Thus, the sum of the diagrams W^; + Wc contains only UV divergence. This IR cancellation 
does not depend on the scheme we use to regularize IR divergences, e.g. by putting a cutoff in 
the tj— integration or by giving a fictitious mass to the gluon. 

Notice, that the IR and UV poles in Wc have opposite coefficients, thus one can formally 
set Wc = 0. In this case however the pole of VF^ for D = A has to be interpreted as an UV 
singularity. 

Summing eqs. ()2.17|) . ()2.18|) and the symmetric contributions we obtain the one- loop expression 
for the unrenormalized Wilson loop ()2.13|) 



47r^/2 



C42zf^{v ■ y - z0)]^-^r(3 - D/2)r{D - 3) (" (J^^ + ^) • (2-20) 
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By subtracting the poles in the MS— scheme we obtain 

W^ll = ^Cf (-L^ + L - ^n'^ , L = ln(z(p - ^0)) + 7e (2.21) 

where 

p={vy)fi={Py)^, 2/=(0,0,y_). (2.22) 

From eqs. ()2.20p and ()2.21|) we can directly see that W{Cs) to one loop depends only on the 
variable p and the possible singularities are in the upper half plane 

W{Cs) = W{p - tO) . (2.23) 

This is due to the fact that p is the only scalar dimensionless variable formed by v, y and /i. 
The "— zO" prescription comes from the position of the pole in the free gluon propagator in the 
coordinate representation. Moreover, expression ()2.2H1 implies that under complex conjugation 

W{p-i{)) = {W{-p-iQ)* . 

It is this property that ensures the reality condition of the distribution function ()2.12|) . 

3. Two- loop calculation 

In this section we perform the two-loop calculation of W{Cs)- In order to simplify the calculation 
and reduce the number of diagrams to compute we use the nonabelian exponentiation theorem 
|14j . According to this theorem we can write 

oo oo 

W^l + Y, (^) VW = exp^ (v)"^^"^ 

n=l n=l 

where w'-"'^ is given by the contributions of Pi^^") with the "maximal nonabelian color and fermion 
factors" to the n— th order of perturbation theory. At one loop we have only the color factor Cp- 
At two-loops the maximal nonabelian color factor is CaCf and the fermionic factor is CpNf. 
Here, Ca is the quadratic Casimir operator in the gluon representation and Nf the number of 
light quarks. From this theorem we have W^-^^ = w^^^ and 

iy(2)=l(^(l))2+^(2) 

and the contributions with the abelian color factor Cp are contained only in the first term. 

In fig. 6 we list all nonvanishing the two-loop diagrams which contain color factors CfCa and 
CpNf. As for the one loop case, the paths in these diagrams do not contain the two rays from y 
to cxD because their contributions cancel due to causality of Wilson lines in ()2.10p . 

The color factor of abelian-like diagrams of figs. 6.1-6.6 is Cf{Cf — ^Ca)- These diagrams 
contribute to w'^'^^ with the color factor —^CfCa- Diagrams with gluon self-energy of figs. 6.7 and 
6.8 will contain contributions with both color factors CfCa and CpNf. The latter one comes from 
the quark loop contribution. Finally, diagrams of figs. 6.9-6.11 are of nonabelian nature involving 
three-gluon coupling and their color factor is proportional to CfCa- We have omitted diagrams 
with abelian color factor Cp, diagrams vanishing due to the antisymmetry of three-gluon vertex, 
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diagrams proportional to y^ = and self-energy like diagrams obtained by iterating the one-loop 
diagram in fig. 4c. The reason for neglecting of the latter type of diagrams is the following. As in 
the one-loop case the self-energy diagrams have both IR and UV poles which cancel each other. 
On the other hand in the sum of all two loop diagrams the IR singularities cancel completely. 
Therefore, as we did for one-loop case, we neglect the self-energy diagrams and interpret the IR 
poles in the diagrams of fig. 6 as UV singularities. 

For the diagrams of figs. 6.2-6.5 and figs. 6.8 and 6.10 we have to consider contributions 
involving cutted propagators. As we have observed before, there is no difference between real 
and virtual gluon propagating in the time-like direction {D{z) = D^{z) for time-like z). This 
gives rise to a partial cancellation between real and virtual contributions. This cancellation has 
been already used to simplify the analysis in the one-loop case ()2.19|1 and we will further exploit 
it in this section. It turns out that to compute singular contributions for D = A we can replace 
the cutted propagators in fig. 6 by full ones. This is due to the fact that one can always deform 
the integration path in such a way that gluons propagate in the time-like direction. 



3.1. Abelian like diagrams 

The general form of the contribution of the abelian like diagrams of figs. 6.1-6.6 is (a = 1, . . . , 6) 



Wa 



9' 



47r^/2 



CF{CF-'^CA)Mp-ti 



S-2D 



UD) 



(3.1) 



where Ia{D) is a function of D and the index a refers to the diagrams in fig. 6 and p is given in 

(E2S1). 

Diagram of fig. 6.1 gives 

/O i-Q rO rl 

dh dt2 dh dt,D{zs-z,)D{z^~Z2) 
oo Jh Jt2 Jo 



with 



D{z, - z,)D{z, - z,) = (^^^^ 



2 r 



ih - ts)H2it2 - —U] 



l-D/2 



where Zi = vti, i = 1, 2, 3 and ^4 = yt4. This diagram contain UV divergence for zi — Z3 ^ 
corresponding to vertex renormalization. The rest of the diagram has the same divergences as 
one loop diagram of fig. 4a. Namely, it has a cusp singularity for Z4^ — Z2 —>■ and light-cone 
coUinear singularity for Z2 ^ 0. Performing the integration we obtain 



h 



T{D/2 - l)r(7 - 3D/2)T{2D - 7) 
6(3-D)(4-D)3 



Diagram of fig. 6.2 gives 



/O ^0 pO pod 

dh dt2 dh dUD{z3-z,)D+{z4-Z2) 
00 Jti Jt2 JO 



with 



D{z,-z,)D^{z,-Z2)=(^^^^ 



(ii - ^3)^(^2 + ^4 + i0){t2 + U- — + 10) 

V- 



l-D/2 
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:i) 




(2) 




(3) 




(4) 










(11) 

Figure 6: Nonvanishing two-loop diagrams containing the "maximally nonabelian" color CaCf 
and the fermionic CpNj factors. Due to the exponentiation theorem, these are the only diagram 
we need to evaluate to compute the two-loop contribution of W{Cs)- The blob denotes the sum 
of gluon, quark and ghost loops. 

where Zi = vti, i = 1,2,3 and Z4 = y — vt^. The analysis of singularities is similar to the previous 
case. We have a UV singularity for zi — z^ ^ and a IR pole originated from one-loop diagram 
of fig. 4b. The result of the integration is 

T{D/2-l)T{7-3D/2)T{2D-7) 



2(3-D)(4-D)2(5-D) 



Diagram of fig. 6.3 gives 



W. 



-g^CpiCp - '^CA){vyy / dt^ 



dty 



dU 



t2 



dUD{z3 - Zi)D+{z4 - Z2) , 



where Zi = vti, Z2 = yt2, Z3 = yt^ and Z4 = y—vt^. Since all singularity in t4— plane lie at the lower 
half-plane we deform the t4— integration path from [0, 00) to (— oo,0]. After this transformation 
Z4 is replaced hj z'^^ = y + vt^ and z'^ — Z2 becomes time-like vector for which cutted and full 
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propagators are the same. We can than replace 

D{z3 - zi)D+{z4 - Z2) => D{z3 - zi)D+{z'4 - Z2) = D{z3 - zi)D{z'^ - Z2) 
which gives 



D{z^-z^)D{z'^-Z2) 



T{D/2-l] 
47r^/2 



tit2{ti-—h){U + —{l-t2)) 

V- V- 



l-D/2 



This diagram has two hght-cone colhnear singularities for Zi ^ Q and 2:4 -^ y. Evaluating the 
integral we obtain 

T\3-D/2)T\D-3) / _ T\5 - D) ' 

^~ (4-D)4 V ~ r(9-2D) 

The diagram of fig. 6.4 gives 



dti / dt2 / dts / rft4D 

00 Jti Jo Jts 



+ (2:3 -2;i)L'+(2;4-2;2) 



where Zi = vti for i = 1,2 and 2;^ = y — vtj for j = 3, 4. As in the previous diagram we can deform 
z^ and Z4 integration paths by replacing Zj = y — vtj with Zj = y + vtj for j = 3, 4 corresponding 
to the replacement 

D^zs - z,)D+{z4 - Z2) ^ D4z', - Zi)D4z', - Z2) = D{z'^ - z^)D{z'^ - Z2) 

and we have 



D(4 - z^)D{z'^ - Z2) 



V{D/2-l] 



ih -t2 + —){U - t2)ih - ti + ^)(t3 - t^] 



l-D/2 



This diagram has a single IR pole for Z4 — Z2 ^ 00 and z^ — Zi ^ 00 simultaneously. Evaluating 
the integral we obtain 



1 T{2D - 7) 
'A-D2T^{D/2~1] 



+ 0{{4-D) 



Diagram of fig. 6.5 gives 



/O pO pi POD 

dh dt2 dh dUD{zs-z,)D+{z4-Z2), 
00 Jti Jo Jo 

where Zi = vti for i = 1,2, z^ = yt^ and z^ = y — vt^. Deforming the z^ integration path we 
replace z^ with z'^ = y + vt^ 



D{z3 - zi)D+{z4 - Z2) =^ D{zs - zi)D+{z'^ - Z2) = D{z3 - zi)D{z'^ - Z2) 



and we have 



D{z3-z^)D{z'^-Z2) 



T{D /2-r 

47r^/2 



ti{h-—h){t4-t2 + —){t4-t2) 

V- V- 



l-D/2 
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It turns out that this diagram has no singularities for D = A. To confirm this we evaluate the 
integral and obtain 



2 



1 f3T{3-D/2)T{2D-7) 



-T {D/2 1)| 3 + (4 _ /))3(3 _ /^) (^ r(3D/2-5) 

T{5-D)T{2D-7) 2T{3 - D/2)T{D - 3) 



T{D-3) T{D/2-l) Jj ^^^^ ^^ 

Diagram of fig. 6.6 gives 

W, = g^CpiCp - \CA){vyf f dt, f dt^ f dt, f dUD{z, - z,)D{z, - z^) 

J —OD J t\ Jo J t-J, 

where Zi = vti ioi i = 1,2, Zj = ytj for j = 3, 4, and we have 
D{z,-z,)Diz,-z,)=(^^^^^J 



2 r -| l-D/2 

tlt2(— t3-tl)(— i4-t2) 



This diagram has two cusps singularities for z^ ~ zi ^ and z^ — Z2 —> and two light-cone 
coUinear singularities for zi —>■ and Z2 -^ 0. Evaluating the integral we obtain 

2,n/o iA f 1 fT{5-D)T{2D-7) T{7 - 3D/2)T{2D - 7) 



h = T^(D/2- 1, , , ,, , , , , , . 

- ^(2 - D)i3 - D/2nD - 3) (r(5 - D) - ^^72^^!) )) ~ 4^^} + ^^^^ ~ ^^"^ ' 
where ({3) is the Reimann function. 

3.2. Self-energy diagrams 

The general form of the contribution of the diagrams with gluon self-energy in figs. 6.7 and 6.8 
is (a = 7, 8) 

Wa = (^) Cj. ((3D - 2)Ca - 2{D - 2)Nf) [2^ip - iQ)f-'^h{D) . 

The one-loop correction to the gluon propagator in the Feynman gauge in the coordinate repre- 
sentation is given by 

which differs from the free propagator in the power of z"^ — iO. For time-like vector z the cutted 
one-loop propagator coincides with the full propagator in ()3.2|) . As in the one loop case this 
allows us to treat the sum of these diagrams with all possible cuts by using the full propagator. 
We obtain 



W7 = g^CF{vy) I dh f dhD^'\z2 

J~oo Jo 



Zl 
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where z\ = vti and Z2 = yt2 and 

dh / dt2D^l\z2-z,) 
■00 Jo 

where Zi = vti and Z2 = y — vt2. By performing the integration we get 

T\D/2-l)T{5-D)T{2D-7) ^ T\D/2 - l)T{5 - D)T{2D - 7) 



I7 



16(4 - L))3(l - L))r(L) - 2) ' ^ S{A-DY(l-D)V{D-2) 



3.3. Diagrams with three-gluon vertices 

For the diagrams of figs. 6.9-6.11 containing three-gluon vertex we have the following general 
expression 



Wa = Ic^Cp (-^\ ■ [2^{p - iQ)r'^UD) 



2 

where a refers to the various diagrams in fig. 6 with three-gluon vertex. We simplify the analysis 
of the diagrams with various cuts by computing only the contributions which are singular for 
D ^ 4. In this case we can treat all gluons as virtual. To show this observe that for the cutted 
diagrams at D = 4 we have only infrared and light-cone collinear singularities. The infrared 
singularities cancel. The collinear singularities appear when gluons propagate along the light- 
cone, i.e. when the intermediate point z^ lies on the segment [0, y]. In this case cutted propagators 
coincide with the full propagators. Notice that cusp singularities appear when all gluons interact 
at small distances. They are present only for diagrams of figs. 6.9 and 6.11 for Zi -^ 0. In this 
case all gluons are virtual. Therefore, in the following we study only the contributions to Wa in 
which all gluons as virtual. 

For the diagram of fig. 6. 9 we have 



W^ = Iq^CaCf I dh [ dt2 f dh fd''z,v^^y^'y^-'r^,,,,,{z,,Z2,z,)l[D 

J —OD Jo Jt2 J j_1 



3 

Zi — Z4^) 



where Zi = vti, Z2 = yt2 and Z3 = yt^. By using the expression for the three-gluon vertex in the 
Appendix A we find 

V'y'-y'-T,,^,J.u^. .,) = .(..) (.|^ - y^^ = .(=*) (|^ - |^) ■ (3.3) 

Due to this particularly simple form of the three-gluon vertex the integration over t2 or t^ becomes 
trivial. For the first term the integration over t2 gives the contributions from the end points Z2 = 
and Z2 = Z3. For the second term the integration over t^ gives the contributions from z^, = Z2 
and Z3 = 1. In all contributions, the integral over the intermediate point Z4 is factorized into the 
following expression 



JyZi, Z2, Z3) 



3 

/rf% 1[D{z,-z,) 
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valid for z^ and 23 on the light-cone {z\ = z\ = {z2 — z^)"^ = 0). The pole at Z^ = 4 corresponds to 
a light-cone singularity as Z4 approaches the segment [0, y]. The meaning of the integral over the 
parameter s is the following. For D = 4 one integrates over the free gluon propagator between 
the points Zi and z = SZ2 -|- (1 — 5)2:3. Since the z lies on the light-cone between Z2 and Z3 
the vector z — Zi is time-like. This confirms the expectation expressed at the beginning of this 
subsection that all coUinear gluons propagate in the time-like direction. 
Performing the remaining integration we obtain 

r(5 - D)r(2D - 7) / T^{D/2-l) AT{7-3D/2)T{D/2-l) T\D/2-l) 



4(4-L))3 \{4:-D)T{D-3) 3{4:-D)T{5-D) T{D-2) 

For the diagram of fig. 6.10 one deforms the integration path over z^ = y — vts into z^ = y + vt^, 
replaces cutted propagators by full ones and obtains 

dt, / dt2 / dh / d^z,vf'^y>''n'^'T^,^,^^{z^,Z2,Zs)l[Diz,-z,), 

■00 Jo Jo J j^;^ 

where Zi = vti, Z2 = yt2 and Z3 = y + vt^. By using the three-gluon vertex we find 
.-,-n-r,,,,,3(^i, Z2. zs) = ^{y^^- y|-) - ^{vy) (.A - .|- 

The first term leads to a contribution which is regular for D = 4. To see this notice that the 
singularities arise when the gluons are propagating along the light-like vector y. However, this 
configuration is suppressed by applying the operator y-^ to the gluon propagator. The expression 

y^D{zi — Zi) is proportional to y{zi — Z4) and vanishes for Zi — z^ parallel to y. The second term 
is similar to the one of the previous diagram. We have contributions from the end-points zi = 
and Z3 = y. For Zi = (2:3 = y) the two vectors zi and Z2 {z2 and z^) lie on the light-cone and 
we can apply the identity ()3.4|) . Performing the remaining integrations we obtain 

r(5 - D)T{2D - 7)T{D/2 - 1) fT{D/2 - 1) r(7 - 3D/2) \ 
J 10 — 



\2D - 7)T{D/2 - 1) fV{D/2 - 1) r(7 - 3L>/2) \ 
2(4 -D)4 V r(D-3) r(5-D) J 



For the last diagram of fig. 6.11 we have 



/O ^0 I'oo /■ J^ 

dh dt2 dh d''ziV^^v^^y^-'V,,^^,^{z,,Z2,z,)\{D[ 
-00 J ti Jo J j_-i 

where Zi = vti, Z2 = vt2 and Z3 = yt^. By using three-gluon vertex we obtain 

f d d \ ( d d 

v^'v^'y^'r^,^,^,{zi,z2,z3) = -Hyw--yTr-]+^ivy)[vjT--v- 



Zi — Z4) 



dzi dz2j \ dzi dz' 

d d \ f d d \ 

-i ( y^ 2/7^ + i{vy) Kr - 7^ ■ (3.5) 



dzi dz2j \dti dt 



(-2 



For this diagram we have the following singularities: a cusp singularity for Zi —>■ with i = 
1, . . . , 4; two independent light-cone collinear singularities for zi, Z2 -^ and Z4 approaching the 
segment [0,y]; an ultraviolet singularity from Z2,Z4 -^ z\. 



The operator y-^ in the first term suppresses propagation along hght-hke vector y of gluon 
from point ^4 to Zi or to Z2. This imphes that the contribution of this part of three-gluon vertex 
contains only a triple pole in 4 — D. The second term in (13. 5|) is similar to the one in (j3.3|) for 
diagram of fig. 6.9. We have two end-point contributions with Z2 = and Zi = Z2- For ^2 = 
the two vectors Z2 and z^ lie on the light-cone and we can apply the identity ()3.4|) . For zi = Z2 
the integral is similar to the one of the gluon self-energy correction. 

We finally obtain 

_ r(5 - D)T{2D - 6) [ V{?,D/2 - 5)1(3 - D/2) 1 ^^^;\ 
^"" 4(4-/}) V 3r(Z}-2)(4-D)3 2^^'^V' 

Recall that all expressions in this subsection are valid up to terms which are regular for D = 4. 

3.4. Renormalization at two-loop order 

Since the diagrams of fig. 6 have nested ultraviolet divergences corresponding to the renormal- 
ization of the vertices and propagators one has to include additional counter-terms. In the 
MS— scheme their contributions are given to two loops by 



w. 



(2) 
c.t. 



/«aV r^V ^N\mn ,n^l4-i.£(3^::^/2)^CD-3) / i i\ 



To obtain the final expression for the renormalized w to two-loops we add the counter-terms, 
subtract the poles in the MS— scheme and take into account combinatorial factors. By using the 
non-abelian exponentiation theorem w^^-' is obtained by omitting the colour factor C'jp in ()3.1|) . 
The final contributions to w'-^^ from the various diagrams (a = 1, . . . , 11) have the following form 



w 



(2) ^ / <-^ ^ 2 



(^) Cf [CA{AaL^ + 5,L3 + CaL^ + DaL) + Nf{E^L'' + F^L^ + G,L) + 0(LO 



where L is given in 1)2.211) and for the various diagrams the nonvanishing coefficients are given 
by 









Ae = 
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_ 1 
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_ 1 

18' 














Bi = 
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~~ ~9 


Br = - 


-I Bg -- 


1 
3' 


Bn = 


- -i B 
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.t. - 


_ 11 
18' 










c,= 


1 
3' 


C2 = l, 


Cs = - 


2 
n 

6 ' 


Ce = - 


72^ ' ^^ 
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31 
36 
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Cs 


_ 5 

6' 


C,- 


_ 25 2 1 
144 " 2 ' 
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Gn = 
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144 " 6 ' 




C,,- 
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11 
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72" 
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3' 
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D4 = h 
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^7 = -S- 
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16 


■v^^ 


Ds 
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31 
36' 


D, = ^C(3) - I - 
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- -K(3), 


Du 


= -m- 


^ + iC(3)- 


1 
6' 




Dct. 
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144" 










Ey = |, 


Ect. - 


1 

9' 




















F7 = 
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_ 1 
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. = -T2^'- 












Summing 


all contribut: 


ions we finally obtain 
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(3.7 
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where 



B = -%C^ + \Nf, 



D = (!C(3)-^vr^-^)C^+(^vr^-^)iV,. 
This expression has the following properties: 

• the coefficient of L^ vanishes; 

• the coefficient of L^ is proportional to the one-loop beta-function; 

• the interpretation of the remaining coefficients in front of L^ and L will be clear from the 
RG equation for light-like Wilson loop we are going to discuss in sect. 4. 

The two-loop calculation con&ms the analytical dependence in fl2.23|) which means that all 
singularities of W{Cs) in the complex p— plane lie in the upper half plane. After integration over 
?/_ in ()2.12|) this leads to the spectral property 

F(x,/i/M) = forx>l. (3.8) 

It has been observed jT^j that the presence in ()3.7p of the L^ term seems to be in contradiction 
with the renormalization properties of Wilson loops. Assuming that W{Cs) is renormalized 
multiplicatively, the L^ term should vanish. However one should notice that we are dealing with 
a light-like Wilson loop which has additional light-cone singularities. Moreover, the fact that the 
coefficient of L^ is just the one-loop beta-function suggests that light-like Wilson loop obeys a 
RG equation. This equation is discussed in the next section. 

4. Renormalization group equation for Wilson loop on 
the light-cone 

To evaluate the structure function for x —>■ 1 given in eq. ()2.12|) one needs to compute W{Cs) in 
()2.13|) to all orders in perturbation theory. A powerful method to resum the expansion is the use 
of renormalization group equation. In this section, following ref . ^j , we deduce the RG equation 
of Wilson loops with path partially lying on the light-cone. 

Recall that for x away from 1 the operator product expansion on the light-cone allows us 
to relate the //—dependence of the structure function with the ultraviolet properties of local 
composite twist-2 operators obtained by expanding in powers of y- the matrix element in ()2.1|) . 
This dependence is described by the evolution equations [S] 

fig- + P{g)-^]F{x,fi/M) = J dzP{x/z)F{z,fi/M), (4.1) 

where the splitting function P{z) is singular for z —>■ 1. For x —>■ 1 the structure function ()2.12|) 
is given in terms of W{Cs), which is a nonlocal operator. This suggests that it is convenient to 
treat the Wilson loop as a nonlocal functional of gauge field rather than to expand it into sum 
of infinitely many local composite operators. This can be done by using RG equation of W{Cs), 
which gives the dependence on the renormalization point fi. Since W{Cs) is a function of the 
single parameter p ~ /i?/_, from its /i dependence we directly obtain the dependence on ?/_ and, 
through Fourier transform in ()2.12|) . the dependence on x for the structure function. 
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4.1. Wilson loop renormalization away from the light-cone 

Renormalization group equation for the Wilson loop away from the light-cone is well known [^ 
and depends on the explicit form of the path. The integration path for W{Cs) is shown in fig. 2. 
It has two cusps at the points and y where the quark is probed by the photon. The important 
property of this path is that the segment £2 lies on the light-cone. As we discussed in the two loop 
calculation, the presence of the L^ term entails that the renormalization properties of W{Cs) are 
different from the ones of Wilson loops with path away from the light cone. 

Suppose for a moment that £2 lies away from the light-cone, i.e. y"^ 7^ 0. Then, the dependence 
on fj, of the renormalized Wilson loop expectation value is described by the RG equation 

L-^ + P(g)-^j lnW^(,2^o)(Cs) = -reusp(7+,^) - reusp(7-,^) , (4.2) 

where 7± are the angles in Minkowski space between vectors ±?/^ and v^ = p^/M with f ^ = 1 

V, ^ (^^) 

cosh 74- = ± — -=^ . 

yy 

The cusp anomalous dimension rcusp(7,5') is gauge invariant function of the cusp angle JH]- In 
the limit of large angles 7 we have [T7] 

rcusp(7,fl') = 7rcusp(5') + C'(7°) , 

where the coefficient Fcuspls') is known to two loop order and is given by 

r..(,)^^C..e)^C.(c.(|-g-.,l). (4.3) 



For the integration path of fig. 2 we have y^ = leading to an infinite 7±, thus eq. ()4.2|l becomes 
meaningless. One can directly check that the two-loop result of previous section does not satisfy 
equation ()4.2|1 . However, as shown in subsect. 4.2, there is a simple way to find the generalization 
of the renormalization group equation for the light-cone Wilson loop. 

4.2. Wilson loop renormalization on the light-cone 

Renormalization group equation for Wilson loops on the light-cone has been proposed in PH] 
and can be obtained as follows. First, one slightly shifts the integration path away from the 
light-cone by setting |/^ 7^ and keeping the cusp angle 7± = |ln(4(t>j/)^/j/^) large. Since 7± 
is a logarithmic function of {vy) and y"^ we need to define its analytical continuation away from 
positive {vy) and y"^. The proper expressions for 7± can be deduced from one- loop calculation of 
Wilson loop ^H] and are given by 

A{{vy)-iQf _ A{{vy) + iQf 



7+ = ^ In 5 ^ > 7- = 5 In 



2 ^2_^o ' '"2 ^2_^o 

where the "— iO" prescription comes from the position of the singularity of the free gluon prop- 
agator in the coordinate representation. By using these expressions we now differentiate the 
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renormalization group equation (j4.2|) with respect to the variable (vy) 



d d\ d 



iiiW(Cs) = -r, 



cusp 




2rcusp(5') 
(vy) 



iO 
(4.4) 
The two terms originate from two cusps of the path of fig. 2 which he on the opposite sides of 
the cut. This is the reason for the appearance of complex conjugation in the second term. 

The variable y"^ disappeared from this equation and one can formally set y"^ = 0. This is the 
proposed renormalization group equation for the Wilson loop on the light-cone ^HI- One easily 
check that two- loop expression (|3.7p of the previous section does satisfy eq. (j4.4p . 

Notice, that Wilson loop on the light-cone depends only on a single variable L = ln(i(p — 
iO)) + 7e. Thus equation ()4.4|) becomes very powerful since one can integrate it and obtain RG 
equation for the light-like Wilson loop 



U^ + f^(9)jA W{Cs) = - [2r,,,p((7)L + r(^)] W{Cs 



(4.5) 



where T{g) is the integration constant. From this equation we can see the origin of the various 
terms in the two-loop calculations in ()3.7|) . The appearance of the one- loop beta function in front 
of the L^ term is obvious from (j4.5p . The coefficient of L^ is proportional to a sum of rcusp(fi') 
and one-loop beta-function. The coefficient of L is given by 
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(4.6) 



While rcusp(fi') is a universal number, T{g) depends on the path under consideration. The unusual 
feature of the RG equation ()4.5p is that the anomalous dimension given by the coefficient of 
W{Cs) in the r.h.s. depends on p, i.e. on the renormalization point /i and ?/_. Actually this 
property leads to the evolution equation for the structure function as we shall discuss in the next 
section. 



5. Evolution equations 

The RG equation for the distribution function near the phase space boundary is obtained by 
using the representation ()2.12|1 and properties ()2.22|1 and ()2.23|1 . Introducing the dimensionless 
parameter a = [Py) we can write 



(5.1) 



/OO J 
— e'"(^-")iy(a/i/M - iO) 
-OO 2vr 

with the inverse transformation 

H{fi/M)W{afi/M -iO)= I dx e~*"(^-^)F(a;, p/M) , 

^— OO 

where the range of x— integration takes into account the spectral property of the structure func- 
tion in (j3.8|) . The renormalization group equation in (j4.5|) gives 

V F{x, fi/M) = L^ + Pig) ^] F{x, fi/M) = j^ dz P{1 + x - z)F{z, fi/M) (5.2) 
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where 

/OO 7 
e*-(i-) {-2reusp(^) ln[z(a/i/M - zO)e^-] - !((?) + V \nH{fi/M)} . 
-OO ^TT 

From the analytical property of the integrand we have 

P{z) = for 2 > 1 , 

this is the reason for setting to x the lower limit for z in ()5.2j) . To compute P{z) we use the 
representation 

io a 
After a careful treatment of the a ^ singularities, we obtain 

'0(1-2 
T 



P(2) = 2r,,3p(^) ( ^— ^ ) +5(l-z)%) (5.3) 



where (...)+ is the standard plus-distribution and 

hig) = -Tig) + 2r,,sp(^) In — + P In H{^/M) . (5.4) 

Due to the factorization theorem 0, P(2;) should not depend on fi and the same is then true 
for the function h{g). This means that the /i dependence in the term involving the coefficient 
function H{fi/M) should be compensated by a contribution from W{Cs)- One can consider now 
fj5.4p as the RG equation for the function Hi^n/M). While only soft gluons contribute to rcusp(5'), 
both soft and collinear gluons (and quarks) contribute to h{g). Therefore the function h{g) is 
not fixed by the RG equation of W{Cs) and should be directly computed [IH] from Feynman 
diagrams. 

Notice that the function F{x,fi/M) defined in ()5.H) satisfies the evolution equation ()5.2|) for 
any value of x. However, this function has the physical meaning of the quark distribution only 
for X ^ 1. 

The evolution equation (j5.2|) in the soft limit a; ~ 1 and 2; ~ 1 can be written in the standard 
form ()4.H) where P{z) is the two- loop quark splitting function for z near 1. From ()5.3|) the general 
form oi P{z) is given by 

{l-z)P{z)=2T,^,p{g)eil-z). (5.5) 

Using the expression ()4.3p for rcusp(fi') one easily verifies that the result of two-loop calculations 
[TSj obeys this relation. 

The evolution equation (|5.2j) for x ~ 1 is obtained from the RG equation for the light-like 
Wilson loop. Because of the universal structure of RG equation, one finds that any distribution, 
which can be represented in terms of the light-like Wilson loops, satisfies equation ()5.2j) with the 
kernel P{z) defined by ()5.5|) . This is the case for the quark and gluon structure and fragmentation 
functions (see subsect. 2.3) at large x. For the gluon distributions one should replace the colour 
factor Cf by Ca in (IO|l . 

Let us consider the moments of the function F{x,fi/M) defined in (j5.1|) . 

Fn{fi/M)= / dxx''-^F{x,fi/M) = H{fi/M) / dxx""'^ / —e'''^^-''^W{an/M-iO). (5.6) 
Jo Jo J -00 27r 

23 



For large n the integral over x receives the leading contribution from the a; — > 1 region in which 
F{x,^/M) coincides with the quark distribution function. By using the two-loop expression 
for W{aiJ,/M — iO) in (|3.7|) we obtain the corresponding Fn{^/M). Before this we note that 
W{an/M—iO) is given to any order in perturbation theory by a sum of powers of L = ln[z(o"/i/M— 
iO)] + je- Hence, in order to find Fn{fi/M) one needs to expand the following basic integral 

/.....-. /;je--K4-.o)]--^(A)-'^^(A„)-',,.,v„,,, 

in powers of 6. This relation implies that for large n and arbitrary 6 the integral is given by the 
integrand evaluated at a = —in. Performing this trick in ()5.6|) we get the simple expression 

F^{fi/M) = H{fi/M)W{-infi/M) (1 + 0{l/n)) , (5.7) 

which means that the large— n behavior of the moments of the quark distribution function is 
given by the Wilson loop expectation value evaluated along the path of fig. 2 with formal iden- 
tification (Py) = —in (or p = —inn/M) in ()2.22|) and ()2.23p . From ()5.7|) and the non-abelian 
exponentiation theorem we obtain that large— n corrections to F„ exponentiate. From the one- 
and two-loop results in (|2.21|) and (|3.7|) we get that the exponent of F„ contains a sum of powers 
of In n up to as In n terms in one loop and a^ In n terms in two loops. The sum of all these 
corrections to all orders can be done by using the RG equation fl4.5|) for light-like Wilson loops. 
From this equation and from ()5.7p we find that the large— ra behaviour of Fn{p/M) is governed 
by 

^IT- + P(9)i-] Fnifx/M) = - [2reusp(^) Hne^^fx/M) + T{g)] F^M) . (5.8) 



dn dg ^ 

To test the relation ()5.7|1 we differentiate Fn{fi/M) w.r.t. /i and substitute eqs. ()4.5j) and ()5.4j) to 
obtain 

^W' + ^^^^^) ^"^^/^) = " [2rcusp(^) ln(ne^-) - h{g)] F^M) . (5.9) 

One easily checks that the same equation follows from ()5.2|) with the l.h.s. equal to the moment 
of the evolution kernel (|5.3p for large n. 

6. Concluding remarks 

To conclude we would like to mention a close relation between the analysis here presented and 
the heavy quark effective field theory (for a review see ref. jTHj) which seems to be a powerful 
tool for analyzing of heavy meson phenomenology. The emission of gluons from heavy quarks can 
be treated by the eikonal approximation. This implies that the propagation of the heavy quarks 
through the cloud of light particles can be described by Wilson lines and that the effective heavy 
quark field theory can actually be formulated [20] in terms of the Wilson lines. Therefore, the 
renormalization properties of Wilson lines discussed in this paper are related to the ones for the 
effective theory. For example one finds that the "velocity dependent anomalous dimension" is 
the cusp anomalous dimension (j4.3p . 

The central point of our analysis was the RG equation (j4.5j) for the generalized Wilson loop 
expectation value, W{C), with path partially lying on the light-cone. The cusp anomalous 
dimension rcusp(5') entering into this equation is a new universal quantity of perturbative QCD 
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which controls the behaviour near the phase space boundary of hard distributions. The evolution 
equation for these distributions corresponds to the RG group equation for W{C), moreover the 
splitting function near the phase space boundary is related to the cusp anomalous dimension. 
This equation does not imply that W{C) is renormalized multiplicatively. 

It is well known that the structure function of deep inelastic scattering gets large perturbative 
corrections for x ^ 1 which need to be summed 0IS|- They come from two subprocesses: from 
quark distribution function and from cross section of the partonic subprocess. In the present 
paper we considered the x ^ 1 behavior of the distribution functions. We were able to control 
in eqs. ()5.8|) and ()5.9p their large perturbative corrections to all orders using the renormalization 
properties of light-like Wilson loops. The same method can be applied for the summation of 
large perturbative corrections to the partonic cross section jS 0]. It will be described in the 
forthcoming paper. 

Note added on February 4, 2005 

Recently the two-loop calculation of the Wilson loop W(Cs) has been performed by E. Gardi in 
ref. |2Tj. The obtained results for the Feynman integrals Ji, . . . , /lo agree with our expressions 
while III differs by ~ 1/{D — 4) term. We repeated the calculation of In and reproduced the 
result of [21]. In the present version, we updated expressions for the integral In, Eq. fl3.6|) . the 
D— coefficient, Eq. ()3.7|) and the anomalous dimension T{g), Eq. ()4.6|) . by taking into account 
the contribution of the additional term. 

We are most grateful to E. Gardi and M. Neubert for useful discussions and correspondence. 

Appendix A: Feynman rules in the coordinate represen- 
tation 

In this appendix we recall the Feynman rules for calculation the generalized Wilson loop expecta- 
tion value in the coordinate representation using dimensional regularization. The D— dimensional 
free gluon propagator in the Feynman gauge is given by D^^{x) = —g^^D{x) where 

^ ' J (27r)^ F + iO 47r^/2 ^ ' 

The "cutted" propagator D'^{x) = —g^^D^{x) associated to a real gluon is defined as 

D+{^) = J £^e-^'^27re{ko)6{e) = ^^^Jj^~^\ -2{x^ - zO){x. - zO)]'-^/' 

where the last equality holds only for x^ = 0. We note that for gluon propagating in the time-like 
direction (x+ > 0, a;_ > and xy = 0) cutted and full propagator coincide. For the three-gluon 
vertex with three gluon propagators attached we have 

^^1^2^3(^15 ^2, Z3) J d Z4^ Y[i=l D{Zi — Zi) 

= -i (^^1^2 (9f - 9f ) + ^'^^W (^m _ 5^1) + ^MiM3 (5M2 _ ^M2)) J ^D^^ |-J3^^ ^(^^ _ ^.) 
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For the gluon attached to the point on the ray ii of the integration path of fig. 2 and propagating 
to z we have igv^j_ dtiD{vti — z) . The analogous expressions for a gluon attached to the 

segment £2 and the ray £3 are igy^ f^ dt2D{yt2 — z) and —igv^ f^ dt^D{y — vt^ — z) respectively. 
The scalar products and integration measure in terms of light-cone variables are 

a± = —^{ao ± 03), a = (ai, 02), (ab) = a^h^ + a+6_ — a ■ b, d a = da+ da^ d ^^a 

V2 

for arbitrary D— dimensional vectors a^ and b^. 
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